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BOUNDARY  PROBLEM  OF  SUPPORTING  SURFACE  IN  NONSTATIONARY  LIMITED 
FLOW 


B.  N.  Belousov,  V.  I.  Orishchenko 
(Irkutsk) 


In  this  article  we  examine  a method  of  solving  the  boundary 
problem  for  a wing  of  arbitrary  aspect  ratio  which  makes  periodic 
arbitrary  oscillations  of  infinitely  small  amplitude  in  a limited 
fluid  flow  having  a velocity  at  infinity  of  Vq . 


Given  is  a brief  derivation  of  the  main  integro-dif ferentlal 
equations  of  the  problems  obtained  in  [1]  and  a method  for  solving 
these  equations  and  obtaining  the  numerical  hydroaerodynamic 
characteristics  of  the  problem  is  developed. 


The  boundary  problem  for  the  potential  of  accelerations  0 is 
formulated  as  [2] 


A0  = o 

0+  -o_  = o 
0 — 0 
^(0)  = o 


g^Sp 

g^L 

x—*oo 

seSo. 


Here  fi  is  the  space  occupied  by  the  fluid;  Sp  - projection  of 
supporting  surface  in  direction  of  undisturbed  flow;  L - trailing 

FTD-ID(RS)I-0124-77  1 


edge  of  supporting  surface;  P - operator  whose  form  Is  determined 
by  the  type  of  limiting  surface  Sq. 


The  solution  to  the  problem  gives  integral  operator  Ay, 

2 

assigned  in  space  L^CS^)  with  values  in  space  C (n). 

The  necessary  properties  of  operator  Ay  are  given  in  [2], 
and  [3].  In  the  studied  problem  operator  Ay  takes  the  form  of: 


where 


— + sign  F 


i] 


ds. 


• r^V(x  — ty  + (y  — ny  + (z  — C)*: 
r, = V{x-  ly  + (y- 1)-  + (2  + ^ + 
f — I— solid  screen 

5ignP»=  I 4- 1 — f ree  surface  of  heavy  fluid 


The  integral  equation  of  the  problem 


AtT^^^Fig),  g^S^ 

by  representing  the  solution  in  the  form  of  3 components 


7 = 7i  + 7,  + Ts 

[2],  [3]  is  reduced  to  the  equation  system: 

^o^A,^l  = ^1  + C g^Sp 
AfoiA^^j  — — ikFi 

g^Sp, 

4 


(1) 

(2) 

(3) 


here:  N 


01 


■■i 


[ ] dT; 


YlCC'(3p) 


YjCC’(Sp) 


regular  solution  related  to  presence  of  break  in 
tangential  velocities  during  transition  of  surface  S 
regular  solution  which  describes  inertial  motions; 
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- singular  solution,  determined  from  equation  (3) 
transferred  to  space  of  velocity  potential  by 
operator  N. 

Equation  (1)  is  not  used,  since  it  contains  an  indeterminant 
constant  C for  the  case  of  a supporting  surface  of  arbitrary 
aspect  ratio. 


Thus,  the  inertia  component  of  the  solution  Y2  is  determined 
from  equation  (2),  which  in  form  coincides  with  the  equations  of 
the  stationary  theory,  while  the  vortex  components  y-j^  and  y^ 
determined  from  equation  (3). 

The  kernel  equations  (2)  and  (3) 

Kix,  5,  y. ,,)  = riZEIEZilMIZEj)! 

(jf— e)(>-ii) 

is  represented  by  the  Leydlo-Sobolev  approximation  for  a wing  of 
arbitrary  aspect  ratio  in  the  form  of  two  singular  kernels  with 
singularities  with  respect  to  (x  - 5)  and  (y  - n),  respectively, 
which  from  the  theoretical  standpoint  is  incorrect,  although  it 
does  lead  to  satisfactory  final  results. 


In  this  study  the  approximation  for  this  kernel  contains  a 
regular  term  in  addition  to  the  singular  components,  and  this 
makes  it  possible  to  preserve  the  main  properties  of  the  kernel 
K(x,  C,  y,  n)  [1]. 


+ sign(jc  — E) 


fci-,-! 


Sin  k 


(y ->,)</*]. 


(4) 


where  . , . 

P(X)^0.  X-O:  P (yJ-o.  X-*oo. 

P(X)-.oo,  X — oo;  ^-*0. 
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The  kernels  of  equations  (2)  and  (3)>  which  consider  the 
limited  nature  of  the  flow,  are  approximated  by  functions  which 
have  the  form  of  the  regular  terms  in  the  kernel  approximation 
K(x,  c,  y,  n) . 


The  system  of  equations  (2)  and  (3)  based  on  approximation 
(4)  is  reduced  to  the  system  of  one-dimensional  Integro-dif ferentlal 
equations  of  the  Prandtl  type 


lltP,(k)  f'  , Pi(k)  f (>1)  ^ , I 

- f 0 0-  ^ J 0,  (v-1)  rf,  I . 

-I  -1 


-I 


Here 


G(y  — -n)  = - 


ji  — i, 


{e  +e-*''*- 


-2  +2«]  * 

— e ^ •'  J sin  A (y  — Yj)  rf/t, 


<«  -2P  ^ 

O,  (y-Y,)  = J 


I -2  Fp/'— "N +2h]  * 

+ « L u;  J isink(y-n)dk. 


Prom  equation  (6)  when  k 0 we  get  the  stationary  theory 
equation . 


The  equation  for  the  inertial  component  of  the  solution  is 
represented  in  the  form  of 


+ l +• 

Mill  J Aj(i,y)0(x-^)d(  + -^j'^7(x,yi)0(y-yi)(fyi-^-i^F,; 

-1  ’ 

0 (X  - ^ - J *- 

_ g-2  IP  (1)  + 2*1  *j  * (jc  _ I)  dk. 


4 


(7) 


As  an  example  of  solving  equations  of  this  type  let  us 
analyze  the  method  of  solving  the  equation  of  the  stationary 
theory 


-1 


-I 


(8) 


Here : 


^)  = j T(5,  ■n)di;  r(y,)^  J t(5.  v)  dl 

-r  -I 

The  solution  y(C>  n)  is  sought  in  the  form  of  a series 

1 (£,,)_  )/l^  K 2 2 ■’„£  V. 

n— 0 /7^o 


Using  the  method  of  moments  with  weight  factor  V 1 —y*  l/i  —jc» , 
equation  (8)  is  reduced  to  a system  of  algebraic  equations 

* k 

« S S ^np'^Bpgj  = Ogn  (9) 


where  coefficients  B 


npqs 


t*’ 

«-0  n-O 

are  represented  in  the  form  of  series 


in  powers  of  functional  parameters 

T*  = /4/i»+  1 - 2A, 

^P(x,-yFw+T-p(^), 

Vp,.,-/4[/.  + ^]’+l-2[*  + £2i.]. 

In  order  to  solve  equation  system  (9)  a program  was  compiled 
on  the  BESM-4  computer.  Calculation  results  in  the  form  of  the 
distribution  of  circulation  T over  the  wing  span  are  shown  in 
Fig.  1. 
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Pig.  1. 

The  solution  to  an  analogous  equation  by  the  Multopp  method, 
which  did  not  consider  the  distribution  of  circulation  over  the 
chord,  was  given  in  the  work  by  I.  P.  Tkachenko  [4].  The  results 
of  the  calculations  agree. 

Equation  (7)  for  inertia  component  Y2  solved  by  the  method 

used  to  solve  equation  (8).  The  solution  of  h)  is  sought  in 

the  form  of 

T($.  s S ApIV- 

/i«»o  n—0 

The  weight  factors  are  used  in  the  same  form.  Calculation  results 
are  shown  in  Fig.  2. 


Fig.  2. 
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